Measurements of the magnetic field induced by a turbulent flow of liquid metal 



in 
o 
o 

(N 
O 

Q 

(N 



C/5 

O 



> ■ 
in ■ 

: 

(N . 
O ■ 

tn : 

o . 

o : 

^ : 

Or 



X 



M. D. Nornberg, E. J. Spence, R. D. Kendrick, C. M. Jacobson, and C. B. Forest)] 

Department of Physics 
University of Wisconsin-Madison 
1150 University Ave. 
Madison, WI 53706 
(Dated: December 23, 2005) 

Initial results from the Madison Dynamo Experiment provide details of the inductive response 
of a turbulent flow of liquid sodium to an applied magnetic field. The magnetic field structure is 
reconstructed from both internal and external measurements. A mean toroidal magnetic field is 
induced by the flow when an axial field is applied, thereby demonstrating the omega effect. Poloidal 
magnetic flux is expelled from the fluid by the poloidal flow. Small-scale magnetic field structures 
are generated by turbulence in the flow. The resulting magnetic power spectrum exhibits a power- 
law scaling consistent with the equipartition of the magnetic field with a turbulent velocity field. 
The magnetic power spectrum has an apparent knee at the resistive dissipation scale. Large-scale 
eddies in the flow cause significant changes to the instantaneous flow profile resulting in intermittent 
bursts of non-axisymmetric magnetic fields, demonstrating that the transition to a dynamo is not 
smooth for a turbulent flow. 

PACS numbers: 47.65.+a, 9L25.Cw 
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I. BACKGROUND 

The generation of magnetic fields by flowing 
electrically-conducting fluids is a long-standing problem 
in plasma physics, astrophysics, and geophysics. Of par- 
ticular interest is the role of turbulence in either support- 
ing, or suppressing, magnetic fields. Dynamos have been 
created in the laboratorMrby driving helical flows of liq- 
uid metal through pipeslra Measurements of the onset of 
magnetic field growth agree with predictions from lami- 
nar theory, suggesting that turbulence played little role 
in these experiments. The pipe geometry, however, limits 
the development of eddies to the width of the pipe, a scale 
much smaller than the device size, thereby inhibiting dy- 
namics due to large-scale turbulence. Since astrophysical 
flows lack this scale separation, understanding the role of 
large-scale turbulence is especially important in model- 
ing the dynamo. The Madison Dynamo Experiment was 
built to characterize these dynamics. 

The experiment, shown in Fig.^j, uses two impellers to 
generate a double- vortex flow of liquid sodium in a 1 m 
diameter spherical vessel. Liquid sodium is used because 
of its high conductivity (a = 10 7 O _1 m _1 ). The choice 
of a spherical geometry is motivated by the computa- 
tional work of Dudley and James, which demonstrated 
that simple time-stationary vortices, can generate mag- 
netic fields at relatively low speeds U The flow has been 
measured in hydrodynamic experiments and is predicted 
by laminae dynamo theory to generate a transient mag- 
netic fieldQ by a slow-dynamo mechanism that stretches 
and twists-.magnetic field lines to regenerate the initial 
seed field.cl Due to the low kinematic viscosity iv — 



7 x 10 -5 m 2 /s), and hence low Prandtl number of liquid 
sodium (Pr = [Lqctv = 9 x 10 -4 ), and the large magnetic 
Reynolds numbers required for spontaneous field growth 
(Rm = noaLVo ~ 10 2 ), the flows required for a dynamo 
are extremely turbulent (Re = LVa/v = Rm/Pr ~ 10 5 ). 
These types of flows are also studied in experiments at 
Maryland] and CadaracheE The open geometry of these 
experiments provides the opportunity to study the role 
of large-scale turbulence in magnetic field generation. 

The purpose of this correspondence is to provide a de- 
scription of the Madison Dynamo Experiment and to re- 
port the initial measurements of the magnetic field in- 
duced by the flowing liquid metal. The theory of laminar 
kinematic dynamos and hydrodynamic measurements of 
the flow are reviewed in Sec.[n|. The experimental appa- 
ratus and its diagnostics are described in Sec. [II. Mea- 
surements of the mean magnetic field induced when an 
external magnetic field is applied to the turbulent flow 
are compared to predictions from the laminar theory in 
Sec.|y|. Measurements of the magnetic field fluctuations 
are used to study the properties of MHD turbulence in 
Sec.0. 



II. THE KINEMATIC DYNAMO 



A. Theory 
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The evolution of the magnetic and velocity fields in 
an incompressible conducting fluid is governed by the 
magnetic induction equation and Navier-Stokes equation 



FIG. 1: (a) Photograph and (b) schematic of the Madison Dynamo Experiment. The sphere is 1 meter in diameter. It is filled 
with 105-110° C liquid sodium, and a flow is created by two counter-rotating impellers. Two sets of coils, one coaxial with and 
one transverse to the drive shafts, are used to apply various magnetic field configurations. The magnetic field induced by the 
flow is measured using Hall-effect sensors both on the surface of the sphere and within tubes that extend into the flow. 



with a Lorentz forcing term: 
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where /io is the vacuum magnetic permeability, a is the 
fluid conductivity, p is the fluid density, and v is the kine- 
matic viscosity. The time scale for magnetic diffusion is 
T cr = PQ<Ja 2 '', where a is a characteristic size of the system. 
Recast in dimensionless units, the induction equation be- 
comes 



dB 

~dt 



RmV x V x B + V 2 B, 



(3) 



where the magnetic Reynolds number Rm = [ioaavo is a 
measure of the rate of advection compared to the rate of 
magnetic diffusion (^o is a characteristic speed) . The rel- 
ative importance of the Lorentz force in the flow dynam- 
ics is given by the interaction parameter N ~ aaB 2 /pvo, 
where Bq is a characteristic field strength. For a kine- 
matic dynamo, N <C 1 so that the velocity field evolves 
independently of the magnetic field. If the flow is sta- 
tionary (i.e., the flow geometry is constant in time), the 
induction equation becomes linear in B; it can be solved 
as an eigenvalue equation by expanding B as 



B(r,<) = ^B,(r)e 



(4) 



where Xi are the growth rates of the magnetic eigenmodes 
Bj(r). A dynamo is produced when at least one eigen- 
value has a positive real growth rate. 

There are two fundamental requirements for a flow to 
produce a dynamo. The first requirement is a sufficiently 
fast flow speed so that the advaetion of the magnetic 
field overcomes ohmic dissipation.13 There is a minimum 
magnetic Reynolds number Rm cvl t below which resistive 
diffusion dominates the evolution of the field. As Rm — > 



i?m cr it, it is expected that the flow will more effectively 
amplify the initial seed field to produce a dynamo. The 
second requirement is feedback — the induced field must 
reinforce the initial seed field. Although the flow may 
induce a strong response from an initial seed field, there 
is no feedback to continue the cycle of magnetic field 
generation if the induced field is perpendicular to the 
seed field. 

The amplification and feedback can be quantified by 
the gain, defined as 



B„ 



gam 



1 COS S + ^applied 



B, 



applied 



(5) 



Here, -Bmduccd is a measure of the amplitude of the in- 
duced field, ^applied is a measure of the amplitude of 
the applied field, and 5 is the angle defining the rela- 
tive orientation of the induced field to the applied field. 
The mechanism which produces a dynamo is desc ribed in 
terms of amplification and feedback in Sec. [I B and the 
use of Eq. |^ in analyzing meas urements from the experi- 
ment is described in Sec. IV B 



B. Experimental flow modeling 

Although the impeller-generated flow in the experi- 
ment is turbulent, the mean flow can be approximated by 
a laminar model. Flows in a spherical geometry are most, 
easily modeled by the Bullard and Gellman formalism.^! 
The velocity field is described by a spherical harmnnic 
expansion of toroidal and poloidal stream functional^ 



V(r, 9, 4>) = J2 (V x V x [ s e , m (r)Y e m (8, 4>) r 



V x 



^W^Wif]) (6) 
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FIG. 2: (color) (a) Velocity field profiles fitted to LDV measurements of water flows generated by impellers identical to those 
used in the sodium experiment. The impeller rotation rate is 16.7 Hz, corresponding to Rmtip = 100 based on impeller tip 
speed. A contour plot of the toroidal flow is shown in the upper hemisphere, and poloidal stream lines are shown in the lower 
hemisphere. The poloidal flow is directed inward at the equator and outward at the poles. The two rectangles represent the 
location of the impellers. The sr,o, S2,o, S4,o, ti,o, t2,o, and £4,0 profiles are used for the flow model, (b) The probability 
distribution function for an LDV measurement of the toroidal velocity in the water model of the sodium experiment. The 
measurement location is shown as the circular dot in (a). 



and the magnetic field is described by an expansion of 
flux functions 



B(r,9i) = ^(VxVx [ S e , m (r)Y e m (6, 0) f 



Lm 



V x 



^ >ro (r)F/"(M) r]) , (7) 



where the Y™(9, <f) are the spherical harmonics, si^ m (r) 
and t£ yTn (r) are radial scalar profiles describing the 
poloidal and toroidal velocity field, and Si -m (r) and 
Ti >m {r) are radial scalar profiles describing the poloidal 
and toroidal magnetic field. The experimental flow is 
axisymmctric and is composed of primarily the £2,0 and 
S2,o terms. Hence, it is called a t2s2 flow. 

The velocity field depicted in Fig.||(a) is determined 
from Laser Doppler Velocimetry (LDV) measurements 
performed in an identical-scale water model of the 
sodium experiment a Liquid sodium at 120°C has the 
same kinematic viscosity and density as water at 40° C, so 
the flow measured in the water model should correspond 
to the sodium flow. The magnitude of velocity fluctua- 
tions is extremely large, as can be seen from the proba- 
bility distribution function from LDV measurements de- 
picted in Fig.||(b). The flow profile depicted in Fig.||(a) 
is therefore only realized in a mean sense. The instan- 
taneous flow profile can deviate significantly from the 
mean, implying that the eigenmode growth rate should 
fluctuate on the timescale of the flow evolution. The in- 
duced field would then exhibit transient behavior, espe- 
cially for flows near Rm cvlt . The auto-correlation time of 



the LDV measurements of the flow is r c = 60 ± 20 ms for 
Rmti v — 100, which serves as an estimate of the timescale 
for changes in the flow profile. 

An eigenmode analysis of the flow profile shown in 
Fig. ||(a) suggests that a magnetic field should be gen- 
erated for Rm > 175 as seen in Fig.||(a). The structure 
of the magnetic eigenmode with the largest growth rate 
for the t2s2 flow is dominated by the S% t i term in the ex- 
pansion of Eq. ^, corresponding to a dipole field oriented 
perpendicular to the axis of symmetry of the flow (note 
the field lines in the first panel of Fig.|j) . The gain, shown 
in Fig.||(b), is calculated from Eq.|[, where Bmduccd is the 
strength of the transverse dipole field induced by the flow 
subjected to a transverse dipole field of strength -B a ppiied- 
The azimuthal angle between the induced and applied 
fields is 5. It can be seen from Fig.|](b) that the gain 
increases as Rm — > i?m cr it, though the feedback is in- 
sufficient to produce a dynamo until Rm > Rm cr i t . It 
should also be noted that Cowling's anti-dynamo theo- 
rem does not apply.-since the eigenmode breaks the sys- 
tem's axisymmetry.El 

The physical mechanism of the gain and feedback for 
the t2s2 flow is demonstrated by following the evolution 
of a field line in the frozen-flux approximation (the dif- 
fusion term in Eq.|| is neglected). In Fig.||, two field 
lines parallel to the equatorial plane are stretched to the 
poles (located where the drive shafts enter the sphere) by 
the poloidal flow and twisted by the toroidal flow back 
into their original direction. The induced field enhances 
the initial seed field, thereby increasing its magnetic flux. 
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FIG. 3: (a) The kinematic growth rate of the dominant eigen- 
mode versus Rm. The growth rate is calculated from the lin- 
earized induction equation using the flow profile constructed 
from LDV measurements in a water model of the sodium ex- 
periment, where Rm t i p = 100 (Rm = 98 based on maximum 
velocity). It is assumed that the impellers generate the same 
flow profile as the impeller tip speed is increased and that the 
maximum flow speed scales linearly with the tip speed, (b) 
The corresponding gain as defined in Eq.|| The amplifica- 
tion of magnetic flux grows as the flow reaches Rm CT i t = 175. 
Above Rm CI it , the gain is undefined since the kinematic model 
does not account for the saturation of the magnetic field. 

This process continues as tension in the field lines builds. 
The frozen-flux assumption eventually becomes invalid; 
the field line tension presumably is relieved through mag- 
netic reconnection due to resistive diffusion. Once the 
magnetic field becomes sufficiently strong, the interac- 
tion parameter becomes large. The Lorentz force creates 
a torque on the flow, referred to as the back reaction, 
which halts the field growth. 

The laminar analysis above may be inadequate for 
describing the experiment since turbulence also induces 
magnetic fields. Eddies in the flow distort the magnetic 
field, creating small-scale magnetic field structures which 
can affect the large-scale induced field. Mean Field The- 
ory suggests that velocity field fluctuations can enhance 
the resistivity of the fluid by—effectively increasing the 
rate of diffusion (the /3-effect),E£l generate large-scale cur- 
rents through helical fluid motion (the a-effect),0 and 
reduce the magnetic field within the flow due to inho- 
mogeneities in_the magnitude of turbulent fluctuations 
(the 7-effect).E3 The usual assumptions of homogeneous, 
isotropic turbulence and scale separation of fluctuations 
from the mean flow are not necessarily satisfied in the 
experiment; hence the effects described are only used for 
conceptual understanding. 

III. DESCRIPTION OF THE EXPERIMENT 

The experiments are performed in a 1 m diameter, 
1.6 cm thick, stainless steel sphere shown in Fig.[l]. Liq- 
uid sodium is transferred pneumatically from a storage 



vessel, located in a vault beneath the floor, to the sphere 
using pressurized argon. The sodium fills the sphere from 
the bottom until the liquid level rises to an expansion 
tank connected to the top of the sphere. The expan- 
sion tank accommodates changes in the volume of the 
sodium due to variation in temperature. The weight of 
the sodium in the storage vessel is monitored during the 
transfer to determine the liquid level in the sphere. Elec- 
trical contact switches in the expansion tank provide a 
redundant means of determining the sodium fill level. 

The conductivity of liquid sodium varies with tempera- 
ture. It has a maximum value of 1.5 x 10 7 (ilm)~ 1 near the 
freezing point of sodium (98°C) and decreases by about 
4% for every 10°C. Since Rm cx a, the sphere is kept at a 
temperature of 105-1 10°C to optimize Rm without freez- 
ing the sodium. The sphere's temperature is maintained 
by a heat exchange system which runs heat-transfer oil 
through a series of copper tubes on the sphere's surface. 
The system provides 12 kW of heating and 75 kW of cool- 
ing to the oil. Kaowool insulation reduces the sphere's 
ambient heat loss. The heat introduced to the sodium 
by the rotating impellers is removed through the surface 
of the sphere by the heat exchange system. 

Two 30.5 cm diameter impellers generate the flow. 
They are driven by 75 kW motors controlled by vari- 
able frequency drives (VFDs). The motors rotate the 
impellers at rates from 3-30 Hz (Rmtip is between 18- 
180). Since i?m cr it depends on flow geometry, the im- 
pellers have been modified to produce the desired flow: 
Kort nozzles limit the radial thrust and fins on the ex-, 
terior of the nozzles increase the mean toroidal flow.tl 
The magnetic Reynolds number Rm t i P is approximated 
using the rotation rate of the impellers, which is mea- 
sured by digital encoders. Although convenient, this es- 
timate is generous in that the speed of the impeller tip 
is much larger than that of the bulk flow. The VFDs 
record torque and power information. The motor power 
follows a cubic relationship with rotation rate as shown in 
Fig. ^. Once a dynamo is achieved, it is expected that the 
back reaction due to the Lorentz force will increase the 
power required to drive the impellers beyond the cubic 
extrapolation shown. 

At atmospheric pressure, the rotating impellers cause 
cavitation by creating a rapid drop in the local pres- 
sure near the impeller blades, forming bubbles in the 
fluid. When the bubbles collapse against the blades, they 
emit ultrasonic noise which is monitored by a transducer 
mounted to the sphere. Since these bubbles disrupt the 
flow through the impellers, the cavitation must be sup- 
pressed. This suppression is accomplished by pressuriz- 
ing the sphere with argon gas. The required pressuriza- 
tion is determined empirically by increasing the sphere 
pressure until the ultrasonic noise is minimized. For ex- 
ample, to operate at a rotation rate of 20 Hz, the sphere 
must be pressurized to 550 kPa (80psi). 

Two sets of coils generate magnetic fields used to study 
the inductive response of the flow. One set is coaxial with 
the drive shafts, the other orthogonal, as shown in Figfl]. 
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FIG. 4: (color) The mechanism of field generation in a laminar double- vortex flow is modeled in the frozen flux limit (Rm — > oo). 
A field line directed through the equator is stretched by the poloidal flow towards the pole and then twisted back onto itself 
by the toroidal flow. 

fluctuations due to eddies down to the resistive dissipa- 
tion scale (see Sec.|y|). The stainless-steel tubes encasing 
the internal sensor arrays vibrate when the impellers are 
driven at rotation rates above 15 Hz (Rm t i p = 90). Since 
the amplitude of the vibrations increases with flow speed, 
experiments are limited to the lower rotation rates to pre- 
vent damaging the tubes and risking a breach. A compar- 
ison of data from experiments before the internal arrays 
were installed with data from experiments with the tubes 
indicates that the disturbance in the flow due to the tubes 
has negligible effect on the large-scale induced magnetic 
field. Future experiments will be performed without the 
internal sensor arrays to reach higher rotation rates. 
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FIG. 5: The motor power follows the characteristic cubic rela- 
tionship with the impeller rotation rate which is proportional 
to Rm. The curve is extrapolated to the maximum motor 
power available, yielding -Rm max = 150. 



A DC power supply provides the coils with 600 A to gen- 
erate fields up to 100 G. The coils can produce dipole, 
quadrupole, transverse dipole, and transverse quadrupole 
field configurations. For flow speeds of lOm/s, the inter- 
action parameter is N = 10~ 3 ; hence the magnetic field 
is advected passively by the flow. 

The magnetic field is measured using Hall-effect probes 
(Analog Devices AD22151 Linear Output Magnetic Field 
Sensors) on integrated circuits with internal temperature 
compensation. The probes saturate at ±170 G. The sig- 
nals are sent through a low-pass filter to reduce the noise 
level. Since the sphere shields the probes from frequen- 
cies higher than the skin-effect frequency, the low-pass 
filter is designed to have a roll-off frequency that matches 
the skin-effect frequency. For the d = 1.6 cm thick stain- 
less steel sphere (r) ~ 7.2 x 10~ 7 S7m), this frequency is 
/skin ~ {^^o^d 2 ) ~ 1 ~ 1kHz. Probes are positioned on 
a grid on the surface of the sphere capable of resolving 
spherical harmonic modes up to a polar order of I = 7 and 
an azimuthal order of m = 5. Linear arrays of probes in 
stainless steel tubes are inserted radially into the sphere 
and are oriented to measure cither axial or toroidal mag- 
netic fields. Data from the magnetic probes is sampled 
by 16-bit digitizers on PC-based data acquisition cards at 
a rate of 1 kHz per channel, which is sufficient to resolve 



IV. MAGNETIC FIELD MEASUREMENTS 
A. Reconstruction of the mean magnetic field 

The induction effects are studied by applying an axial 
magnetic field to both flowing and stationary sodium and 
comparing the measured fields. The flow's differential ro- 
tation wraps the field lines around the drive shaft axis to 
produce a toroidal magnetic field through the so-called u- 
effect, as seen in Fig.g(a). The contour plot is generated 
by fitting the coefficients of the harmonic expansion in 
Eq.J^to toroidal field measurements (up to order £ = 3). 
The w-effect is very efficient at amplifying the applied 
magnetic field; measurements of the toroidal magnetic 
field shown in Fig. ^(b) indicate that the induced toroidal 
field increases with i?m t i p and is larger than the applied 
field for Rmtip = 120. In contrast, Fig.|g shows that the 
axial magnetic field in the equatorial plane is reduced by 
half. This reduction of poloidal magnetic flux can be ex- 
plained by the effect of flux expulsion due to the strong 
poloidal circulation in the £2s2 flow.H 

The mean field structure shown in Fig.^(a) is obtained 
by an averaging procedure which calculates the peak of 
the histogram of magnetic field measurements from each 
probe. This technique is necessary since some of the sig- 
nals have non-Gaussian statistics (discussed in Sec.^). 
The resulting mean field is predominantly due to odd 
harmonics in accordance with the selection rules govern- 
ing the interaction terms of the induction equation.EI 

The magnetic field outside the sphere is reconstructed 
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FIG. 6: (color) (a) A plot of the toroidal magnetic field produced by the w-effect when a 60 G axial field is applied to a t2s2 flow 
with Rm = 100. The contours are calculated by fitting the expected 2i, m (r) profiles to measurements of the internal toroidal 
field. The dots indicate the position of the Hall probes. External field lines are shown based on a vacuum-field expansion fit to 
measurements from the surface probe array, (b) Toroidal magnetic field measurements near the toroidal maximum above the 
impeller and (c) axial magnetic field measurements in the equatorial plane for various Rm. 
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FIG. 7: (color) Time series plots of the axisymmetric modes 
evaluated at 9 = 0. A constant external dipole field is applied 
and the motors are turned on at t = 7s. 



from the field harmonics calculated from the surface ar- 
ray measurements. The poloidal flux lines are shown in 
Fig. ||(a). Since the Hall probes on the surface of the 
sphere lie outside regions containing currents, the mea- 
sured field can be described in terms of a vacuum field 
potential which satisfies V 2 <£> = 0, where B = — V<f>. 
Solving Laplace's equation for the potential in a spheri- 
cal geometry yields 

*(r, 6, cp) = J2 (ct, m r* + 5pr) YT(6, 4). (8) 

The Ci m terms are due to currents in the external field 
coils whereas the Di m terms are due to currents in the 
flow. The field can be separated into applied and induced 



fields whose radial components are 

B applicd -r = -Y^C^-HYPiO,®, (9) 

B induccd -r = ^^(|+l)y™(0,0). (10) 

The external measured magnetic field can be completely 
described in terms of these expansion coefficients. To 
determine the coefficients Di >m for a given set of magnetic 
field measurements, the applied field is subtracted from 
the measured field, and a design matrix Ay is constructed 
which satisfies B.; = AyDj, where Bj = B r (n,0i,(j>i) is 
an array of measurements of the induced field and D 7 = 
Dtj.mj is the array of expansion coefficients in Eq. (|1C|). 
An example of an element of the design matrix is 

The coefficients are obtained Jsv matrix inversion using 
Singular Value Decomposition^ to solve Dj = (Ay) _1 Bi. 
Axisymmetric mode time series are shown in Fig.[?], 
where the modes have been evaluated at 9 = 0. The 
dipole component of the induced magnetic field grows 
non-linearly with Rrritip and cannot be explained by in- 
duction effects from an axisymmetric flow. The dipole 
response results from a turbulent EMF generated by cor- 
related turbulent fluctuations in the flow.til 



B. Measuring gain with a transverse applied field 

The proximity of the flow to magnetic self-excitation is 
studied by applying a magnetic field perpendicular to the 
drive shaft axis. Recall from Sec.O that the anticipated 
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FIG. 8: (a) Amplification of a magnetic field applied perpen- 
dicular to the drive shaft axis as a function of Rmtip. Gain 
is defined in Eq.|^. The solid line is the gain predicted by the 
kinematic model shown in Fig.|(b). (b) Azimuthal angle be- 
tween the induced transverse dipole field and the applied field. 
Error bars on both plots indicate RMS fluctuation levels. 



structure of the magnetic field generated by the dynamo 
is a dipole oriented in this direction. Figure ||(b) sug- 
gests that such a field should be amplified by the flow. 
As in Sec.||, the gain in Eq. || is calculated using the 
i = 1, to = 1 expansion term in Eq.^ for -Bmduccd, and 
S is the azimuthal angle between the applied and induced 
fields. The gain is less than one if a component of the 
induced field is anti-parallel to the applied field, indicat- 
ing attenuation, and is greater than one if the fields are 
in phase. The gain is determined from measurements 
of the induced magnetic field and is shown to increase 
with i?TOtip in Fig.H(a), though not as quickly as antici- 
pated from the kinematic model. Figure ^(b) shows that 
the induced field is somewhat out of phase with the ap- 
plied field at low i?TOti p and that the alignment improves 



shown in Figg(c) from a probe near one of the drive 
shafts, however, has normal statistics. Examination of 
the structure of the magnetic field during the bursts re- 
veals a strong m = 1 mode shown in Fig. Ifj. The turbu- 
lent flow thus induces an axisymmetric field on average 
punctuated by intermittent bursts that break the sym- 
metry. 

Eddies at the smaller scales create structure in the. 
magnetic field down to the viscous dissipation scale.liS 
This structure is evident in the spatial spectrum in 
Fig. [n]. The spatial spectrum is constructed from the 
power spectrum of the signal from a Hall probe that 
measures the toroidal field just above one of the im- 
pellers. The £;~ 5 / 3 scaling can be derived from the in- 
duction equation using a weak-field approximation. The 
magnetic field is described in terms of a mean and a weak 
fluctuating part, B = Bo + B. The Fourier transform of 
Eq.|3] becomes 



ix + — I B k , w = «k x v k ,w x B , (12) 
Moo"/ 



where second-order advection terms have been assumed 
to be negligible due to the weak-field approximation. 
Magnetic fluctuations are due primarily to the aclvsction 
of eddies by the mean flow (the Taylor hypothesisE-3) and 
so the dispersion relation is approximately u(k) ~ kvo. 
Dividing Eq. n3 by this dispersion relation gives 



iu> k 
kv /iQcruo 



Bl 



. k Vk )U 

l T x — 

k v 



X Bn 



(13) 



At scales k <C k a = ^.quvq = Rm/a, the dissipation term 
is small. Thus, magnetic fluctuations at these scales are 
primarily due to advection of the mean field by velocity 
fluctuations, giving B% u /Bq ~ v t^/ v o- For isotropic, 
homogeneous turbulence, the velocity spectrum is the 



as -Rmtip increases. These measurements will be used Kolmogorov spectrum E K (k) 



Jk (x fc~ 5 / 3 . The re- 



to obtain a flow profile in the sodium experiment which 
maximizes the gain. 



V. PROPERTIES OF MHD TURBULENCE 

The preceding analysis documents the large-scale in- 
duced mean field. Since the flow is turbulent with 
Rm 1 and the magnetic field is relatively weak, eddies 
can twist and stretch field lines at many scales. While 
smaller eddies are advected by the mean flow, larger ed- 
dies can distort the large-scale flow itself, and cause inter- 
mittent behavior in the induced fieldEj Figure [| shows 
signals from two Hall probes on the surface of the sphere 
and their probability distribution functions (PDFs). In 
addition to random fluctuations with a Gaussian dis- 
tribution, the Hall probe signal in Fig.^(a) has large- 
amplitude bursts which skew its PDF shown in Fig.^(b). 
The time average of the signal for this probe is there- 
fore not the most-probable value of the field. The signal 



suiting magnetic spectrum is Eu{k) = B^^/k ot k~ 5 ^ 3 
as seen in Fig.[ll] for the range fco < k < k„. This sc&L 
ing was observed in the Maryland dynamo experiment ,t3 
but in the von Karman sodium (VKS) experiment, a k~ l 
scaling was observed.0 The discrepancy was attributed to 
a saturation of the induction mechanism. 

For the range k 3> k a , the advection and diffusion 



terms become comparable. From Eq. 13, we have the 
scaling B\ u jB^ ~ (fiQd/k) 2 v\ u . Hence, the magnetic 

spectrum is Em{^) oc k~ 2 EK(k) oc fc~ n / 3 in the resistive 
dissipation range, as seen in Fig.^ for k > k a . The 
^-11/3 sca li n g was observed in the VKS experiment, but 
the Maryland experiment observed a steeper spectrum 
due to shielding effects of the stainless-steel vessel. 

The dissipation scale k„ is evident from the knee in 
the wave number spectrum of Fig.|ll|. The results are 
summarized in Tab.|, showing that k a increases with 
i?TOti p . Consequently, the magnetic field gains struc- 
ture at smaller scales as Rm t i p increases, down to scale 
sizes of t a — 2ir/k a — 16 cm at Rm, t i p — 100. A magnetic 
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FIG. 9: (a) Measurement from a single Hall probe on the surface of the sphere near the equator, (b) The probability distribution 
for the signal in (a). The impeller rotation rate for this example is 16.7 Hz (Rm = 100) and a 60 G dipole field is applied. A 
Gaussian fit to the right side of the distribution is shown (dotted line) to illustrate the asymmetry, (c) The time series from a 
probe near the drive shaft axis (or pole) and (d) its probability distribution is shown for comparison. 




-35 



-11 flrfd[Gf 11 



35 



FIG. 10: (color) The mode structure of the intermittent 
bursts corresponds to a dipole field aligned perpendicular to 
the drive shaft axis, as seen from the surface magnetic field 
(the axis of rotation is vertical) . The mean magnetic field has 
been removed to isolate the structure of the fluctuation. 



Reynolds number for the turbulent part of the flow can 
be constructed assuming that k a — Rmt ur bko- Using the 
values in Tab.||, it can be shown that i?mt U rb = 0.2i?mti p . 
This relation implies that the small-scale turbulent eddies 
produced by the mean flow are slower than the injection- 




FIG. 11: (color) The inferred spatial spectrum constructed 
from the Fourier transform of the signal from a Hall probe in- 
side the sphere for three different values of Rm t i p . The probe 
is located just above the impeller near the toroidal maximum 
in Fig.|^. Fluctuations are assumed to be due to convection 
of spatial variations in the field. The dispersion relation is 
lj = kvo where vo is determined from velocity measurements 
in the water model of the sodium experiment. The wavenum- 
ber corresponding to the scale size of the experiment is labeled 
as ko . The wavenumber corresponding to the dissipation scale 
ka is labeled for the three spectra. 
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TABLE I: The resistive dissipation scale determined from the 
spatial spectrum in Fig.[ll]. The velocity is measured in the 
water model of the experiment at the probe's position. 



Rmup 


«o [m/s] 


K [m- 1 ] 


30 


0.39 


20 


60 


0.78 


29 


100 


1.2 


40 



scale eddies. The dissipation scale is therefore only about 
one decade smaller than the injection scale rather than 
the two decades expected for a flow with Rm ~ 100. As 
a result, the small-scale dynamics of the flow has little 
role in generating magnetic fields; it is the largest eddies 
in the flow which induce magnetic fields and which have 
the greatest impact on the transition to a dynamo. 



VI. CONCLUSION 

The magnetic field induced by a turbulent flow of liq- 
uid sodium demonstrates flux expulsion and the cj-cffcct. 



The flow generated in the experiment induces a field that 
provides amplification and feedback, the necessary in- 
gredients of a dynamo. Power spectrum measurements 
show that the resistive dissipation scale is inversely pro- 
portional to Rm. The turbulent flow creates intermittent 
bursts of large-scale magnetic fields that have a structure 
similar to the magnetic field predicted from the laminar 
kinematic dynamo model. Consequently, as a flow nears 
-Rwcriti the transition to a dynamo appears to be inter- 
mittent rather than smooth. 
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